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Abstract 

MnW04 is regarded as a canonical example of multiferroic materials, where the multiferroic ac- 
tivity is caused by a spin-spiral alignment. We argue that, in reality, MnW04 has two sources of the 
spin-spirality, which conflict with each other. One is the Dzyaloshinskii-Moriya (DM) interactions, 
reflecting the P2/c symmetry of the lattice. The P2/c structure of MnW04 has an inversion cen- 
ter, that connects two Mn sublattices. Therefore, from the viewpoint of DM interactions, different 
Mn sublattices are expected to have opposite spin chirality. Another source is competing isotropic 
exchange interactions, which tend to form a spin-spiral texture with the same chirality in both 
magnetic sublattices. Thus, there is a conflict between DM and isotropic exchange interactions, 
which makes these two sublattices inequivalent and, therefore, breaks the inversion symmetry. Our 
theoretical analysis is based on the low-energy model, derived from the first-principles electronic 
structure calculations. 
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I. INTRODUCTION 



The idea of breaking the inversion symmetry by some complex magnetic order has at- 
tracted a great deal of attention. It gives rise to the phenomenon of improper multiferroicity, 
where the ferroelectric (FE) polarization is induced by the magnetic order and, therefore, 
can be controlled by the magnetic field. Alternatively, one can change the magnetic structure 
by applying the electric field. This makes multiferroic materials promising for the creation 
of the new generation of electronic devices.-*^ 

One of magnetic textures, that breaks the inversion symmetry, is the spin spiral. More- 
over, this property is universal, in a sense that it can take place in all types of magnetic com- 
pounds, irrespectively of their symmetry, as long as the the spin-spiral order is established.- 
Probably, this is the reason why the idea of the spin-spiral alignment became so popular 
in the field of multiferroics: many properties of such compounds are interpreted from the 
viewpoint of spin-spiral order," and the search for new multiferroic materials is frequently 
conducted around those with the spin-spiral texture.- 

Nevertheless, the spin-spiral alignment alone is insufficient for having a FE activity: 
although the spin-spiral order breaks the inversion symmetry, one can always find some 
appropriate uniform rotation of spins, which, after combining with the spatial inversion, 
will transform the spin-spiral texture to itself (see Appendix |A]). The FE polarization in 
such a situation will be equal to zero. In order to make it finite (and in order to fully break 
the inversion symmetry), one should ban the possibility of uniform rotations in the system 
of spins. This can be done by the relativistic spin-orbit (SO) interaction, and today it 
is commonly accepted in the theory of multiferroic materials that the spin-spiral alignment 
should be always supplemented with the SO coupling.-i^ However, as long as the SO coupling 
is involved, the spin-spiral texture becomes deformed. Then, we face the question whether 
the FE activity is caused by the spin-spiral order itself or by the deviation from it.- The 
answer to this question is very important for understanding the origin of multiferroicity. 

MnW04 is an important material in the field of multiferroics. Although the transition 
temperature to the multiferroic phase is low and the FE polarization is weak, the material 
is fundamentally important, because: (i) its multiferroic phase has a spin-spiral texture (or 
is believed to be the spin spiral); (ii) the relative directions of the spin-spiral propagation 
vector, the spin rotation axis, and the FE polarization seem to be consistent with predictions 
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based on the spin-spiral theory.®"^" Thus, MnW04 is regarded as one of successful manifes- 
tations of the spin-spiral theory and, actually, the discovery of the FE activity in MnW04 
was guided by this theory. 

More specifically, MnW04 crystallizes in the monoclinic P2/c structure. It exhibits three 
successive antiferromagnetic (AFM) transitions: at Tn3~ 13.5 K, Tn2~ 12.5 K, and Tni~ 7.6 
K. The API phase, realized below Tni, is nearly collinear and the directions of spins alternate 
along the monoclinic b axis as ttii (the reason why the AFl phase is also called as the 
'ttii phase'). The AF2 phase (Tni< TkT^) is incommensurate and is typically ascribed 
to the spin spiral with the propagation vector qAF2 = (—0.214,1/2,0.457) (measured in 
units of reciprocal lattice translations). The spin moments lie in the plane formed by the 
easy magnetization direction in the plane ac and the b axis. The AF3 phase {T]<s2<T<T^s) 
is characterized by a collinear sinusoidally modulated AFM order with the same q. The 
AF2 phase is ferroelectric. The polarization vector is parallel to the monoclinic b axis: 
P = (0,P{,,0). Pb is relatively small: it takes a maximal value (~ 50 fiC/m^) at around 
Tni, and then monotonously decreases and vanishes at around T^2-^ Other phases display 
no sign of FE activity. 

In this work, we will argue that the multiferroicity in MnW04 has a more complex 
origin. Namely, we will show that: (i) It is not sufficient to have a simple spin-spiral 
texture; (ii) The spin spiral should be deformed by some conflicting interactions, existing in 
the system. In MnW04, these are Dzyaloshinskii-Moriya (DM) interactions and isotropic 
exchange interactions. The conflict of these two interactions breaks the inversion symmetry 
and gives rise to the FE activity in the AF2 phase. 



II. METHOD 

We follow the same strategy as in our previous publications, devoted to multiferroic 
manganites.^ Our basic idea is to construct a realistic low-energy (Hubbard-type) model 
for the Mn 3d bands, which would be sufficient for describing basic magnetic properties 
of MnW04 at a semi- quantitative level. The model is constructed in the basis of Wannier 
orbitals, starting from the electronic structure calculations in the local density approximation 
(LDA) with and without the relativistic SO interaction. The details can be found in the 



review article (Ref. Ill[). All calculations have been performed using experimental parameters 



of the crystal structure.— Since heavy W atoms is an important source of relativistic effects, 
the SO couphng was included already at the level of ordinary LDA calculations, and thus 
obtained electronic structure was used as the starting point for the construction of the low- 
energy model. Thus, although the W sites were not explicitly included to the model, the 
relativistic effects, associated with these sites, were allowed to contribute to the dispersion 
of the Mn 3d bands and, in this way, were taken into consideration during the analysis of 
the low-energy model. All obtained parameters of the model Hamiltonian are collected in 
After the construction, the model was solved in the mean-field Hartee-Fock (HF) 



Ref. 
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approximation. 

We would also like to briefly comment on merits and demerits of our technique in com- 
parison with the first-principles approach: 

(i) It is certainly true that the construction of the model Hamiltonian is based on some 
additional approximations,— and is sometimes regarded as a step back in comparison with 
first-principles calculations. On the other hand, it would not be right to think that the first- 
principles calculations for the transition-metal oxides are free of any approximations. The 
necessity to treat the problem of on-site Coulomb correlations, which is frequently formulated 
in terms of the LDA+U approach (with some phenomenological correction "-|-?7" , borrowed 
from the Hubbard model) make it similar to the model approach: it also relies on additional 
approximations, although of a different type. The typical approximations are the choice of 
the parameter U and the form of the double-counting term, which are both ill-defined in 
LDA-I-L'^.— Thus, at the present stage, it would not be right to say that one techniques is 
more superior than other: rather, they provide a complementary information for the analysis 
of material properties of the transition-metal oxides. As we will see below, results of our 
model analysis are well consistent with available first-principles calculations for MnW04.— 

(ii) The search of the ground state for noncoUinear magnetic textures can be very time 
consuming, even at the level of mean-field HF approximation: the rotations of magnetic 
moments towards new equilibrium positions after including the SO interactions can be very 
slow and require tens of thousands of iterations, as in the case of MnW04. In such a 
situation, the model Hamiltonian approach is very useful, because it allows us to find a 
fully self-consistent magnetic structure, which is not always accessible for the ordinary first- 
principles calculations.- 

(iii) The model Hamiltonian approach can be very useful for the analysis and interpre- 



tation of results of complex electronic structure calculations and the experimental data. 
Particularly, in this work it will help us to elucidate the microscopic origin of the magnetic 
inversion symmetry breaking in MnW04. We will also show that the experimental behavior 
of long-range magnetic interactions in MnW04 (Ref. llTI ) can be naturally understood in the 
framework of the superexchange (SE) theory and reflects similar behavior of the transfer 
integrals, derived for the low-energy model. 



III. RESULTS AND DISCUSSIONS 

First, we solve the low-energy model in the HF approximation, by assuming the collinear 
ferromagnetic (FM) alignment of spins, and drive parameters of isotropic exchange inter- 
actions. For these purposes we employ the theory of infinitesimal spin rotations.— The 
procedure corresponds to the local mapping of the change of the one-electron energy onto 
the spin Hamiltonian of the form l-is = ~ Jij^i ' ^j, where is the direction of spin at 
the site i and the summation run over inequivalent pairs of sites. The results of these cal- 
culations are explained in Fig. [T] Alternatively, one can apply the theory of superexchange 
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FIG. 1. (Color online) (a) Lattice of Mn sites with the notations of isotropic exchange interactions 
(numbered in the increasing order of interatomic distances). Note that the P2/c structure of 
MnW04 has two Mn sublattices, which are shown by different colors and denoted as '/' and '//', 
respectively. These sublattices are transformed to each other by the inversion operation, (b) 
Distance-dependence of isotropic exchange interactions: results of calculations, using the theory 
of infinitesimal spin rotations near the ferromagnetic state (denoted as 'inf') and the theory of 
superexchange interactions with Aex = 5 eV (denoted as 'SE'), in comparison with the experimental 
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(SE) interactions, by considering the energy gain caused by virtual hoppings (t^"^') in the 
second order of perturbation theory— For the configuration of Mn^"^, this expression is 
extremely simple: Jij = — X^mm' ^ii^'^ii'^^/^ex, where Aex is the intraatomic splitting be- 
tween the majority- and minority-spin states and the summation runs over all 3(i orbitals: 
m {m') = xy, yz, 3z^—r^, zx, and x'^—y^. Thus, all SE interactions are expected to be 
antiferromagnetic. Then, by taking Agx = 5 eV and using the values of t^'^' obtained for 



the low-energy model (Ref. Il3l ). one can find that the SE theory excellently reproduces 
the results of the more general theory of infinitesimal spin rotations. Moreover, the value of 
Acx = 5 eV is well consistent with typical estimates in the framework of the Hubbard model, 
where Aex ~ U +4 Jh and the parameters of averaged on-site Coulomb repulsion U and the 
exchange interaction Jh, derived for the low-energy model, are 1.8 and 0.8 eV, respectively.— 
Thus, all estimates are consistent with each other. Therefore, we may conclude that: (i) the 
physically relevant mechanism, responsible for isotropic exchange interactions in MnW04, 
is the superexchange; and (ii) the distance-dependence of Jij refiects the behavior of the 
transfer integrals t^™'. 

Our calculations capture main details of the experimental distance-dependence of Jij 
(Fig. H]).— This is a very important finding, which means that in our model we should be 
able to reproduce the correct magnetic structures of MnW04. The negative aspect is that 
our parameters Jij are 'too antiferromagnetic' (AFM). For example, using the values of Jij 
reported in Fig. [H the Curie- Weiss temperature can be estimated as 6'cw~ Jij/Sk^. It 
yields 6cw= —265 K, which exceeds the experimental value —78 K (Ref. |8|) by factor three. 
The discrepancy cannot be resolved simply by changing the value of the on-site Coulomb 
repulsion U, which is frequently regarded as an ill-defined parameter. Particularly, in order 
to explain the existence of FM interactions, which were observed experimentally in some of 
the bonds, we need an additional mechanism on the top of our model. Such a mechanism 
can be related to the magnetic polarization of the oxygen band, similar to orthorhombic 
manganites.— The first-principles GGA+f/ calculations (where GGA stands for the gener- 
alized gradient approximation - an extension of LDA) also substantially overestimate \Jij\ 
and l^cwl-"^ Although these calculations take into account the effect of the oxygen band, 
the disagreement is probably caused by additional approximations for the double-counting 
energy, which is ill-defined in GGA+U .— Finally, the experimental parameters Jij them- 
selves may be sensitive to the details of fitting of the spin- wave dispersion. In this respect. 



we would like to mention that earlier experimental data were interpreted in terms of rather 
different set of parameters Jij-— Thus, at the present stage, taking into account the com- 
plexity of the problem of interatomic magnetic interactions in MnW04, the agreement with 
the experimental data in Fig. [T] can be regarded as satisfactory. 

Then, we search for the magnetic ground state of MnW04 without relativistic SO cou- 
pling. For these purposes we employ the spin-spiral formalism, which is based on the gener- 
alized Bloch theorem.- The idea is to combine the lattice translations with rotations in the 
spin subspace. It allows us to treat the incommensurate magnetic textures with an arbitrary 
vector q = (g^, 1/2, gc) in the reciprocal space such as if the magnetic primitive cell was of 
the same size as the crystallographic one.'^ The results of these calculations, which were also 
performed in the HF approximation for the low-energy model, are summarized in Fig. |2l 
For q = (—1/4,1/2,1/4), the collinear ttii ^-nd noncollinear spin-spiral states are nearly 
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FIG. 2. (Color online) Results of Hartree-Fock calculations for the spin-spiral states with q = 
(ga, 1/2, Qc): total energies versus qa for different values of qc (in units of 2-k/c). The energy of the 
collinear tti-l state is marked by the symbol (8>. 

degenerate. Nevertheless, the total energy of the spin spiral continues to decrease with the 
increase of g^. Therefore, the 'nonrelativistic' ground state of MnW04 is expected to be the 
spin spiral with qNR ~ (—0.175, 1/2,0.45), which is lower in energy than the ttii state by 
about 0.5 meV/Mn. Moreover, none of these states is ferroelectric: although the spin spiral 
breaks the inversion operation, the latter can be always combined with an appropriate spin 
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rotation, which will transform the inverted spin spiral to the original one. Thus, the FE 
polarization will vanish (see Appendix IXj) . 

Therefore, we should find the answer to the following questions: (i) Which interaction 
stabilizes the AFl structure and make it the ground state of MnW04 at low T? (ii) What 
is the mechanism of the inversion-symmetry breaking in the AF2 state, which yields finite 
Pfc? For these purposes we consider the relativistic SO interaction. It gives rise to such 
important ingredients as the single-ion anisotropy and DM interactions. 

The effect of single-ion anisotropy on the magnetic texture of MnW04 can be studied 
by enforcing the atomic limit and setting all transfer integrals equal to zero.— Generally, 
the single-ion anisotropy is expected to be small for the nearly spherical configuration of 
Mn^+. Since the symmetry operation {C^|c/2} of the P2/c space group transforms each 
Mn site to itself (here, stands for the 180° rotation around the b axis and c/2 is the 
additional translation), the single- ion anisotropy will tend to align the spins either parallel 
to the h axis (| |6) or in the ac plane (g ac). In the latter case, the spins are canted off the a 
axis by about 41°, which is close to the experimental value of 35°.- The easy magnetization 
direction corresponds to the in-plane configuration. However, the energy difference between 
hard (||6) and easy (g ac) magnetization directions is small (less that 0.05 meV/Mn), so 
that in the spin-spiral texture it can be overcome by the energy gain caused by the isotropic 
exchange interactions. 

The P2/c structure of MnW04 has two Mn-sublattices (denoted as '/' and '//' in Fig. [1]), 
which are transformed to each other by the inversion operation. Therefore, the DM inter- 
actions will obey the following symmetry rules:— 

(i) all interactions between different sublattices are equal to zero; 

(ii) the DM interactions in the sublattices 'J' and '/J' are related by the identity d;^^ = 
d^Jji = —dji, for any translation R, connecting the sites within one Mn sublattice 
(note that d is the axial vector). 

From now on, it is more convenient to use the extended notations and recall that each site 
of the lattice i is specified by its position r in the primitive cell and the translation R. 

Moreover, due to the symmetry operation {C^\c/2}, all DM vectors lie in the ac plane, 
i.e., similar to the easy magnetization direction, obtained from the analysis of the single- 
ion anisotropy. The parameters of DM interactions can be evaluated by considering the 
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mixed type of the perturbation theory expansion with respect to the SO couphng and 
infinitesimal rotations of spins near the FM state with subsequent mapping of obtained 
results for the change of the one-electron energy onto the spin Hamiltonian of the form 
S(Kii') R'l'^R ^ ^i?']) constructed for each magnetic sublattice.^^ It yields the following 
parameters of DM interactions (in meV): = (-0.01,0,0.01), = (-0.02,0,0.02), and 
^a+c — (0.01,0, —0.01). Interactions in other bonds are considerably smaller. 

The equilibrium magnetic structure, obtained after switching on the relativistic SO inter- 
action in the coUinear AFl state, is explained in Fig. |3l Due to the single-ion anisotropy, all 
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FIG. 3. (Color online) (a) The magnetic structure AFl, as obtained in the mean-field HF calcula- 
tions for the low-energy model with relativistic SO interaction, and (b) its projection onto the ab 
plane. 

spins are confined mainly in the ac plane and canted off the a axis by about 46° (which is 
close to 41°, obtained in the atomic limit). The DM interactions yield an additional canting 
out of the ac plane. The corresponding rotational force f , acting at some Mn spin from its 
neighboring sites in the same sublattice, is given by f = x ^r]- Since d n = —dji, in 

order to contribute to f , the directions of neighboring spins should also satisfy the condition 
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e n = —eji. In the ttii texture, such a situation takes place for R = a and a+c, but not 
for R = c. Moreover, since both dn and eji lie in the ac plane, the force f is parallel to 
the b axis. Finally, due to the alternation of e^^ in the ttii texture, the directions of such 
forces will also alternate. All in all, this explains fine details of the magnetic structure in 
Fig. [Hb- Very importantly, the isotropic exchange interactions, single-ion anisotropy, and 
DM interactions in the AFl phase do not conflict with each other in a sense that all of 
them, even jointly, continue to respect the crystallographic P2/c symmetry. Therefore, the 
inversion symmetry is preserved and the polarization is equal to zero, even after including 
the SO interaction. 

The effect of SO interaction on the spins-spiral alignment can be best understood for q = 
(—1/4, 1/2, 1/2) (see Fig. H]). In this structure, the easy magnetization direction (e G ac) 



FIG. 4. (Color online) (a) The deformed spin-spiral texture, as obtained in the mean-field HF 
calculations for the low-energy model with relativistic SO interaction for q = (—1/4, 1/2, 1/2), and 
(b) its projection onto the ac plane. 

at certain Mn-site alternates with the hard magnetization direction (e||b) at its neighboring 
sites. Thus, there is a conflict between isotropic exchange interactions and the single-ion 
anisotropy. However, since the latter is small, this magnetic structure can be realized at a 
low energy cost. Then, due to the DM interactions, each site with e||b will experience a 
force from its neighboring sites with e G ac. Since d G ac, this force will be parallel to 
the b axis. Alternatively, each site with e G ac will experience the rotational force from its 
neighboring sites with e||6. Since d G ac, this force will also lie in the ac plane. 
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Furthermore, since = — d^, the forces in different sublattices will act in the opposite 
directions. This means that, if some DM vector d^ in the first sub lattice will tend to support 
the right-handed spiral, its counterpart in the second sublattice (d^) will support the left- 
handed spiral. However, the spin-spiral texture, that is formed by the isotropic exchange 
interactions Jij, is always either right- or left-handed in both magnetic sublattices. Therefore, 
we have a conflict between isotropic and DM interactions: if in one magnetic sublattice, the 
effect of these two interactions is added, in other magnetic sublattice, it will be subtracted. 
All these tendencies are clearly seen in fine details of the magnetic structure, obtained in 
the HF calculations with SO coupling (Fig. HJ. First, the canting of spins, lying in the ac 
plane, appears to be different for the sublattices / and //: 21° and 34°, respectively, relative 
to the a axis. Second, as discussed above, the spins with e||6 will experience the force f, 
which is also parallel to b. However, since d;^ = — dj{, the vectors e and f at the same 
site will be either parallel or antiparallel, depending on the magnetic sublattice. Thus, the 
magnitude of spin magnetic moment with e||6 will also depend on the magnetic sublattice. 
Of course, this change of the magnetic moments is small (about 0.01 %, according to our 
HF calculations). Nevertheless, it does take place and is just another manifestation of the 
conflict between isotropic and DM interactions. 

Thus, the conflict between isotropic and DM interactions in the AF2 phase makes the 
sublattices / and II inequivalent. Since in the P2/c structure, these two sublattices were 
connected by the inversion operation, their inequivalency means that the inversion symmetry 
is broken and the actual symmetry of the AF2 phase is lower than P2/c. Furthermore, the 
homogeneous spin-spiral alignment itself is deformed by the DM interactions and, strictly 
speaking, the AF2 phase is no longer the spin spiral. The magnitude of this deformation 
can be seen in the inset of Fig. HI 

Then, we discuss relative stability of different magnetic phases with the SO interaction. 
In this case, the generalized Bloch theorem is not applicable and there only possibility is 
to work with the supercell geometry.- Therefore, we were able to consider only solutions 
with q = (ga, 1/2, 1/2) and the rational qa= —1/3, —1/4, and —1/5. The results can be 
summarized as follows. Amongst AF2 states, the one with qa= —1/5 has the lowest energy. 
The state with qa= —1/4 is higher in energy by about 0.2 meV/Mn. This behavior is similar 
to calculations without SO coupling (see Fig.|2]). Then, the AF2 state with qa= —1/5 appears 
to be nearly degenerate with the AFl state (the energy difference is about 0.07 meV/Mn, but 
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the AF2 state is still low in energy). Although these numerical values are probably on the 
verge of accuracy of our model analysis, this tendency clearly shows that the SO interaction 
additionally stabilizes AFl state relative to the AF2 one. This seems to be reasonable: the 
conflict of DM and isotropic exchange interactions, acting in the opposite directions in one 
of the Mn sublattices, will penalize the energy of the AF2 state. On the other hand, in the 
AFl state, there is no such conflict and the DM interactions will additionally minimize the 
energy of the ttii structure, obtained without SO interaction. Apparently, the main reason 
why the energy of the AF2 phase is still slightly lower than that of the AFl phase is related 
to the fact that the isotropic exchange interactions are overestimated in our approach (see 
Fig. [T]) and, therefore, the effect of the SO coupling and DM interactions on the magnetic 
structure is underestimated. 

Finally, we evaluate the value of the FE polarization in the AF2 phase, using the Berry 
phase formalism,— which was adopted for the low-energy model. The vector of polarization 
is parallel to the b axis, in agreement with the experiment.- Then, the value of Pb can be 
estimated as 2.0, 3.8, and 4.4 /iC/m^ for qa= —1/3, —1/4, and —1/5, respectively. Thus, 
Pi) is small, but it is well consistent with the fact that the DM interactions are also small 
and can be regarded as a small perturbation to the homogeneous spin-spiral texture, formed 
by isotropic exchange interactions. Nevertheless, we would like to emphasize that the DM 
interactions are essential for deforming the spin spiral, breaking the inversion symmetry, and 
producing the finite value of Pb. 

The experimental polarization is an order of magnitude larger: Pb ~ 50yuC/m^.^ How- 
ever, it can be still regarded as a 'small value' in comparison with many other multiferroic 
systems.- Therefore, it is quite consistent with our main idea that Pb is a result of small 
perturbation of the magnetic structure, caused by the relativistic effects. 

There may be several reason why the experimental value of Pb is larger than the theoretical 
one: 

(i) The experimental Pb may also include some lattice effects, in response to lowering of 
the magnetic symmetry;— 

(ii) A unusual aspect of MnW04 is that the inversion symmetry is broken in the "high- 
temperature" phase, while the low-temperature phase remains centrosymmetric (for com- 
parison, the situation in perovskite manganites is exactly the opposite).- Thus, in addition 
to the FE polarization, one should find some mechanism, which would stabilize the phase 
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AF2 over the phase AFl at finite T. This mechanism may involve some temperature effects. 

(iii) The low-energy model is designed for the semi- quantitative analysis. It is not always 
possible to expect a good quantitative agreement with the experimental data, because some 
ingredients can be missing in the model. In principle, the GGA+U calculations are also 
of a semi- quantitative level, because the value of Pb depends on the adjustable parameter 
U.— In our model analysis, the small value of Pb may be again related to the overestimation 
of isotropic exchange interactions (see Fig. [1]). Then, the deformation of the spin-spiral 
texture, caused by the relativistic effects, is underestimated. Therefore, Pb is also underes- 
timated. From this point of view, the behavior of isotropic exchange interactions, the total 
energies, and the FE polarization is consistent with each other, and the main efforts towards 
quantitative description of the FE polarization and the phase diagram of MnW04 should 
be concentrated on the correct description of isotropic exchange interactions. 

IV. SUMMARY 

In summary, we have provided the microscopic explanation for the origin of FE activity in 
MnW04. The multiferroicity in this compound is caused by the conflict of DM and isotropic 
exchange interactions in the AF2 state. Thus, MnW04 is multiferroic not simply because 
of the spiral structure. It is essential to have conflicting interactions, which deform the spin 
spiral, break the inversion symmetry and, thus, give rise to the FE activity. 

Acknowledgements. This work is partly supported by the grant of the Ministry of Edu- 
cation and Science of Russia N 14.A18.21.0889. 

Appendix A: Nonexistence of ferroelectricity in homogeneous spin-spiral state with- 
out spin-orbit interaction 

In this appendix, we will show that, although the spatial inversion is not the symmetry 
operation of the homogeneous spin-spiral state, it can be always combined with an appro- 
priate rotation spins, which transforms the inverted spin structure to the original one. The 
ferroelectric polarization in such a situation will be equal to zero. The prove is extremely 
simple, but flrst we would like to illustrate the basic idea on a cartoon picture for the one- 
dimensional spin spiral (see Fig. |5]). The basic property of the spin spiral is such that: if 
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FIG. 5. (Color online) Cartoon picture, explaining how the inverted spin spiral can be transformed 
to the original one by the uniform rotation of spins. The inversion center is marked by cross. 
The directions of magnetic moments in the original spin spiral (e) are shown by the dark (blue) 
vectors, and those in the inverted spin spiral (Xe) are shown by the light (cyan) vectors. The 
inverted spin structure can be transformed to the original one by 180° rotation of spins around the 
axis n||(e+Xe), which has the same direction at all sites of the lattice. 

e is the direction of spin at certain magnetic site and Xe is the direction of spin at the 
same site after the inversion of the lattice, one can always define the axis n||(e+Xe), which 
will have the same direction at all magnetic sites. Then, the inverted spin structure can be 
transformed to the original one by the uniform 180° rotation of spins around this axis. 

Now, we will provide a more general prove of this statement. To be specific, consider 
the situation, realized in MnW04 and assume that there are two sites in the primitive cell, 
which are located as r and — r, respectively. Thus, these two cites can be transformed to 
each other by the inversion operation X and the inversion center is located in the origin. 
Then, the translations R will specify the location of all other inversion centers, while the 
atomic position will be given by the vectors J2±t. 

In the spin-spiral texture, the directions of magnetic moments are given by 



^ cos q- (Kir) ^ 
sin q ■ [R ± r) 

V » / 



(Al) 



For our purposes, it is sufficient to consider the situation where all spins rotate in the ah 
plane. The generalization for an arbitrary orientation of the rotation plane is straightforward 
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and will be discussed at the end of this appendix. 

Consider the inversion around an arbitrary center Rq. It transforms an arbitrarily 
taken site R±t to Ro—{R±t—Ro) = 2Rq—R^t. Then, the direction of spin at the 
site 2Rq—R^t will change from e2RQ-R^T to eR±T- Thus, our goal is to find a transfor- 
mation 7^, which would rotate Gr±t back to G2Rq r^t- This transformation is the 180° 
rotation around the axis n||(e2i?„-fi=pT+eH±T)- The corresponding transformation matrix is 
given by 

^ cos(2q-i2o) sin(2q-i?o) 0^ 
n = sin(2q • Ro) - cos(2q ■ Rq) • 

[ " " -1/ 

Since TZ does not depend on R, it corresponds to the uniform rotation of all spins. Without 
SO coupling, this is a transformation to the equivalent magnetic state. Thus, although 
X is formally broken, TZX is the symmetry operation, which transforms the spin spiral to 
itself. Since the uniform rotation of spins R does not affect the polarization, the existence 
of the symmetry operation TZI means that the ferroelectric activity in the homogeneous 
spins-spiral state is forbidden. 

This is a general property of the spin spiral: although we have considered a specific 
lattice geometry, which is more relevant to MnW04, absolutely the same argument can be 
repeated, for example, for orthorhombic manganites, where the magnetic sites are located 
in the centers of inversion. 

For an arbitrary orientation of the spin rotation plane, gr^t in Eq. flAll) should be 
replaced by VeR^T, where T) is just another uniform rotation of spins, specifying the orien- 
tation of this plane. Then, TZ should be replaced ViZV^^ . After that all arguments can be 
repeated again. 
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